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Squeezed correlations of particle-antiparticle pairs were predicted to exist if the hadron masses 
were modified in the hot and dense medium formed in high energy heavy ion collisions. Although 
well-established theoretically, they have not yet been observed experimentally. We suggest here a 
clear method to search for such signal, by analyzing the squeezed correlation functions in terms of 
measurable quantities. We illustrate this suggestion for simulated <j)<j) pairs at RHIC energies. 
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INTRODUCTION 

In the late 90's, it was shown [l| that the hadronic 
mass-modification in hot and dense media could lead 
to a a novel type of correlation between bosons and 
their antiparticles. These squeezed back-to-back cor- 
relations (BBC) are the result of a quantum mechani- 
cal transformation relating in-medium quasi-particles to 
two-mode squeezed states of their free, observable coun- 
terparts [H . This is achieved by means of a Bogolioubov- 
Valatin (BV) transformation linking the asymptotic cre- 
ation (annihilation) operators, (Sk), of the observed 
bosons with momentum = (a;fe,k), to the operators, 
(^k), corresponding to thermalized quasi-particles in 
the medium. This transformation is given by — 
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b^kS-k , being Ck = 
cosh(/fe) and Sk = sinh(/fc). For conciseness, we keep 
here the short- hand notation introduced in Ref. [H , where 
(— fc) denotes an opposite sign in the spacial compo- 
nents of the momenta. Since the BV transformation be- 
tween these operators is equivalent to a squeezing oper- 
ation, the coefficient of this transformation, fki.kji^) = 



ilog 



is called squeezing parameter; w^.. 



m and = -|-ki are, respectively the disper- 

sion relation in terms of the asymptotic mass, m, and in 
terms of the in- medium modified mass, m.^,. A complete 
description of the phenomenon should include a param- 
eterization for m.^, depending on the particles' momenta 
and on their coordinates in the hot and dense system. 
For proposing the means to search for it experimentally, 
however, it suffices to asume a linear relation between 
the masses, i.e., = m ± 5m, as was also considered in 

After the publication in Ref. pj], a similar BBC be- 
tween fermion-antifermion pairs was demonstrated to ex- 
ist , if the masses of these particles were modified in- 
medium. Both the fermionic (fBBC) and the bosonic 
(bBBC) back-to-back squeezed correlations are described 
by analogous formalisms, being both positive correlations 



with unlimited intensity. This behavior is in contrast 
to what is observed in femtoscopy, or Hanbury-Brown 
and Twiss effect (HBT), a quantum statistical correla- 
tions among identical particles. In the HBT realm, two 
bosons with similar momenta are positively correlated, 
with intensity ranging from 1 to 2, whereas two fermions 
with similar momenta are anti-correlated, with intensity 
between and 1. Besides, the fBBC and the bBBC con- 
stitute a direct probe of hadronic mass shift in the hot 
and dense media, contradicting naive expectations that 
such information would vanish at the freeze-out surface, 
leaving no trace on correlations. 

In the remainder of this paper, we will focus on the 
bosonic case only, more especifically, on bosons that are 
their own antiparticles, such as 
particle correlation function is written as 



C2(ki,k2) 



A^2(kl,k2) 
7Vi(ki)iVi(k2)' 



(1) 



where the numerator is the two-particle joint distribution 
and the denominator is the product of the two single- 
inclusive distributions. 

The numerator in Eq. ([Ij is proportional to the expec- 
tation value of the four-operator, (al^^al^^aii^ait^-^) . After 
applying a generalization of Wick's theorem to locally 
equilibrated systems 0, |3| , the complete two-particle dis- 
tribution in such cases can be written as iV2(ki,k2) — 

Wki^ks (a^iaiL^akaaki) = Wki^k^ {al^ak^} {al^ak^) + 

(aki«k2)(ak2«ki) + (4i«L)(«k2aki) • 

For estimating the above expectation values the 
asymptotic operators (a, a^) are first written in terms 
of the ones in- medium, {b,w). These last two opeartors 
diagonalize the full, in-medium Hamiltonian, i.e., Hm = 
Ho + Hi ~ J d^k ilkbl.bk', the asymptotic Hamiltonian is 
Hq = J d^k ujka\ak and Hi is proportional do the mass- 
shift. These thermal averages are calculated by means of 
the density matrix operator, p, as (O) = Tr[pO), where 
p = |exp( - J d^k nkblbk), Z - Trip) and T 

is the temperature. The resulting correlation function, 
given by the ratio in Eq. ([1]), can be written as 
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C2(ki,k2) = 1+- 



|G'c(ki,k2 



|Gs(ki,k2 



Go(ki,ki)Gc(k2,k2) Gc(ki,ki)Ge(k2,k2)' 

(2) 

where the denominators represent the product of the two 
spectral distributions, Gc(ki,ki) = Ni(ki) — i^ki^j^ = 

cjki (O't. O'ki)- In the second term, the numerator is writ- 
ten in terms of the chaotic amphtude, Gc(ki,k2) = 
^Wki J Oka ) , and is originated in the indinstin- 

guibihty of the two identical mesons (either cfxp or 
tt'^tt'', in the current discussion), reflecting their quan- 
tum statistics. In normal conditions the third term, 
whose numerator is proportional to the squeezed am- 
plitude, G's(ki,k2) — ^Wki (o-ki ) I does not con- 
tribute. However, if the interactions in the hot and dense 
medium lead to mass-modification, the third term trig- 
gers this striking particle-antiparticle correlation. The 
three terms in Eq. ([5]) contribute together in the 
case of neutral bosons that are their own antiparti- 
cles, such as (fxp or tt^tt^. For charged bosons , such 
as TT^ or K^, these terms split in two separate corre- 
lation functions, envolving different pairs of particles. 
Thus, the first two terms in Eq.® correspond to the 

HBT correlation, C^{k,M) = 1 + o^WG^^u^y 
for identical particle pairs (tt^tt* or K^K^). On 
the other hand, the sum of the first and the last 
terms leads to the squeezed particle-antiparticle corre- 
lation, C,(ki,k2) = 1 + GAS£}G'(i2,i^2) ^ fo'" particle- 
antiparticle pairs (tt^tt^ or K^K^). 



RESULTS ON SQUEEZED CORRELATIONS 



Initial studies of the problem were performed for a 
static, infinite medium P, Later, it was extended 
to the case of finite-size systems expanding with mod- 
erate radial fiow Q. For the sake of simplicity a non- 



relativistic treatment with flow-independent squeezing 
parameter was considered there, which allowed to obtain 
analytical expressions for both the squeezed and the fem- 
toscopic correlation functions Jd\. However, those studies 
focussed on the behavior of the maximum of the squeezed 
correlation function, Cs(k, — k, m*), in terms of modi- 
fied mass, TO*, for particle-antiparticlepairs with exactly 
back-to-back momenta, ki=— k2=k [5|, |6|. In studies of 
the HBT effect, this investigation would correspond to 
focusing on the behavior of the A-parameter, i.e., the in- 
tercept of the correlation function, for identical particles 
with exactly identical momenta. 

Although important for theoretically understanding 
the finite size and flow effects on the squeezed corre- 
lation function, the systematic study of Cs(k, — k, to*) 
does not represent a practical tool to look for the BBC's 
experimentally. In reality, the momenta of the two de- 
tected particle are never exactly back-to-back and the 
in-medium shift in the hadronic mass is not a quantity 
measurable in the detector. For an empirical search of 
the BBC signal, and considering the non-relativistic con- 
text of Ref. we suggest the following. First, select the 
particle and the antiparticle from the same event, with 
momenta (ki, k2), and combine them to form the pair av- 
erage and relative momenta, respectively, K= ^(ki -f k2) 
and q= (ki — k2). Then, analyze the squeezed correla- 
tion function in terms of these variables, similarly to what 
is done in HBT. The maximual value of the BBC effect 
is reached for exactly back-to-back pairs, ki = — k2 = k, 
being located around Kj^^ 0. Therefore, the squeezed 
correlation function should then be investigated by vary- 
ing Ki2 in the region where it is small, for several values 
of qi2, i.e., Cs(ki,k2) Cs(2ifi2, qi2). 

The squeezed correlation function for (fx/) pairs is ob- 
tained by inserting, in Eq.Q, the squeezing amplitude 
and the spectra extracted from the results in Ref. [5| , and 
rewritten in terms K12 and qi2, respectively, as 
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where K12 ± ^qi2 = kii with i = 1,2 are the individ- 
ual momenta, and considering the region where the con- 
tribution of the middle term (HBT part) in Eq. ^ is 
negligible. The flow-modified radius and temperature in 
Eq.© and dU, are given, respectively, as i?* — R^/TjTl 



and = T- 



as in Refs. [5|,@. We adopt natural 
units, h = c = \, througout this article. 

Extending the analogy with usual procedures adopted 
in HBT, we could think that background could be chosen 
experimentally by combining particle-antiparticle pairs 
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from different events. This choice corresponds to con- 
sider an uncorrelated pair, free from identical particle 
exchange effects. However, this is different in the BBC 
case, since the squeezing factor appears in the denomi- 
nator of the correlation function as well. Therefore, on 
searching for the effect of squeezing, we should consider 
the product of the spectra of the particle and the an- 
tiparticle in the denominator of Cs{2Ki2,qi2), each one 
writen as in Eq. 

Naturally, the analysis in terms of the variable 2K 
would not be suited for a genuine relativistic treatment. 
In this case, a relativistic four-momentum variable can be 
constructed as Q^^^^. = (wi— a;2, ki-|-k2) — (g",2K), first 
introduced in [6]. In fact, it is preferable to redefine this 
variable as Qlf,^ = -{QbackY = 4(^1^2 - K'^Ki_,), since 
its non-relativistic limit is Qli,^ — > (2K)^, recovering the 
average momentum of the pair introduced above. 

In Ref.Q, we presented some introductory results on 
(jxp squeezed correlations in terms of |Kj2 1 and |, but 
most of the plots shown there were obtained by attribut- 
ing precise values to the variables in Eq. ([2]), dS]), andQ. 
In the current analysis, however, we follow the above 
procedure in a more realistic estimate, in which the mo- 
menta are generated in a simulation and then combined 
to form IKj^jI and |qj^2l- The binning in such a simu- 
lation should refiect the finite experimental resolution in 
momentum, and experimental acceptance cuts could also 
be introduced, whenever available. It should be stressed 
that the choice of (f) mesons considered here was made 
as a means to illustrate the proposed method, for two 
main reasons: 0's are their own antiparticles and their 
large mass validates the non-relativistic approximation 
considered here. 
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FIG. 1: (Color online) Generated transverse mass distribu- 
tions (arbitrary normalization) are compared with PHENIX 
data for three sets of fit parameters: i) the freeze-out tem- 
perature Tfo = 109 MeV and flow velocity Pt = 0.77, ii) 
Tfo = 157 MeV and Pt = 0.4, used in PHENIX simulation, 
and Tfo = 140 MeV and /3t =< u >= 0.5, adopted here. 



In the simulation we generated the momenta of the 



mesons with the following prescription. For roughly mim- 
icking the experimental cuts, based on PHENIX data on 
(/>'s [rj, we introduced some simple geometrical selection 
criteria in their generation. This was done by merely con- 
sidering the cuts in azimuthal angle and in the pseudo- 
rapidity, as well as by selecting the momentum region. 
In Fig. [1] we show the transverse momenta generated in 
the simulation, as compared to experimental data points, 
for three sets of parameters corresponding to the temper- 
ature, T, and to the radial flow velocity, (u). All three 
sets are in reasonable agreement with the measured dis- 
tributions. In the remainder of this work, we fix T = 140 
MeV. 

We then combine the momenta of the (fxj) particle- 
antiparticle pairs to calculate the squeezed correlation 
function, as explained above. The top two plots corre- 
spond to the simulated Cs(rn»,qi2), keeping the average 
momentum of the pair in a small interval, |Kj^2 1 ^ 1 MeV. 
In the top most plot no cuts were considered, but in the 
middle plot of Fig. [5] (a) , a rough version of the cuts from 
Ref. 7] were introduced in the simulation. No sensitiv- 
ity to those cuts is apparent. At last, for cross-checking 
the simulation code, we compare the squeezed correlation 
functions estimated with the generated pairs, with that 
obtained by considering pairs with exactly back-to-back 
momenta. This is seen in the plot at the bottom of Fig. 
[2l[a), showing the histogram obtained by attributing ex- 
act values to q^^, and also fixing K^^ = 0, similarly to 
what was shown in Ref. 0, 0, IJ ■ This plot in the bot- 
tom shows close similarity with the first two, as would 
be required. 

Since the strength of the squeezed correlation is ex- 
pected to be significant in the low-jKj^jl region, we also 
suggest to plot Cs(K^2 7Qi2) a function 2K^2: ^oi en- 
larging the average momentum region where the corre- 
lation intensity can be significantly above unity. Fig. 
[2] (b) shows the result for the squeezed correlation func- 
tion from simulation, obtained considering the static case 
((w) = 0), on top, and the case with radial flow {{u) = 
0.5), in the middle. In both, we considered that the parti- 
cles were emitted in a finite interval. At = 2 fm/c, during 
which the emission decreases due to a Lorentzian distri- 
bution in time, \Fs{At)f = [l + {uji+iU2fAtY^ flSSi, 
which multiplies the third term in Eq. ([2|) . This factor re- 
duces the signal by almost three orders of magnitude, as 
compared to a instantaneous emission {At — 0). This can 
be seen by comparing the plot in the middle of Fig. [UJb), 
with the one in the bottom. We note that, although the 
reduction of the strength caused by a finite emission pe- 
riod is dramatic, the intensity of the (f>(j> BBC correlation 
is still sizable, suggesting that its experimental search 
is indeed promising. Another particular emission time 
distribution will be discussed below. 

From Figs. [Sfb) we see that, in the absence of fiow, the 
squeezed correlation grows faster from smaller to higher 
values |qj^2 1 than in the presence of flow. However, this 
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FIG. 2: (Color online) Part (a) shows the squeezed correlation function in terms of the shifted mass, m,, and the pair relative 
momentum, jqijl- Results from simulation (top two plots) are compared with those from calculation, with fixed 7^12 — 
(bottom). Part (b) shows the advertised plots for the squeezed correlation function to be searched for experimentally, with 
m, = IGeV. The striking effect of finite emission times is shown by comparing the middle and the bottom plots in part (b). 
Part (c) shows the inverse width of the BBC function reflecting the radius of the squeezing region, for 7? = 7 fm (top) and 
i? = 3 fm (bottom). 



last one is stronger even in the low jq^^l region, show- 
ing that the presence of flow could enhance the signal's 
intensity over a wider region of the (|Kj2|,|qj2 |)-plane. 
The size of the squeezing region is reflected in the in- 
verse width of the curves as a function of 2|Kj2|, being 
narrower (broader) for larger (smaller) radii. In Figs. 

we compare the squeezed correlation functions con- 
sidering that the system radius is i? = 7 fm, also used 
in all the above calculation, to the case of a system with 
smaller extent of the squeezing region, with R = 3 fm. 

Returning to the discussion about the time emission 
distribution, we should emphasize that we do not know a 
priori how this emission should proceed, mainly in the 
simple model adopted here. On the other hand, the 
lack of knowledge of such distribution, naturally does 
not invalidate the experimental search of the hadronic 
squeezed correlations. The squeezing is a fundamental 
phenomenon already detected in quantum optics and it 
should be empirically observed in relativistic heavy ion 
collisions, if the hadron masses are modified in-medium 
by some mechanism. Therefore, as a common practice 
in Physics, we should look for the effect experimentally 
and, once it is discovered, we could try to explain its time 
emission process by means of a suitable model. Never- 
theless, motivated by an analysis made by the PHENIX 
Collaboration 3 we investigated in Ref. [l^l a dif- 
ferent emission distribution in time, by considering the 
effects of a a symmetric, a-stable Levy distribution, i.e., 
|F(At)p = exp{—[At{uji+uj2)]°'}, on the squeezed corre- 
lation function of K^K^ pairs. This functional form had 
been fitted to two- and three-particle Bose-Einstein cor- 
relation functions. The values fitted to date different val- 
ues of the distribution index, a = 1.0 or a = 1.35, were 



fitted to data, depending on the region investigated of 
the particles' transverse momentum or transverse mass. 
Briefiy summarizing those effects, we concluded that, for 
a = 1, such a time factor acting on the squeezing corre- 
lation function reduces its intensity even more dramat- 
ically than the Lorentzian factor discussed here. How- 
ever, it would still lead to measurable quantities, mainly 
if the emission lasted a short period of time, of about 
At — 1 fm/c. However, if Nature favors the higher value, 
a = 1.35, this would result in a very small deviation from 
unity, not detectable by the method proposed here. Nat- 
urally, in the case of (j) pairs, even for a — 1 and At — 1 
fm/c the intensity would be negligibly small for a Levy- 
type distribution, due to their large asymptotic mass. In 
spite of that, as discussed above, we do not know a priori 
the preferred form chosen by Nature for the emission pro- 
cess, which in itself does not invalidate the experimental 
search of this phenomenon. For this reason and for con- 
tinuing the illustration of the proposed method to search 
for the hadronic squeezed states, in the remainder of this 
paper we attain our discussion to the Lorentzian time 
distribution, comparing it to the instantaneous emission 
process only. 



EFFECTS OF SQUEEZING ON <j><j) HBT 
CORRELATIONS 

We next discuss how the HBT correlation function 
could be affected by the in- medium mass-shift. Contra- 
dicting early expectations that the thermalization would 
wash out any trace of mass-shift in this type of correla- 
tion, we find its in-print in HBT, reflecting the presence 
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of the squeezing factor, fij{m,mf,), in the chaotic am- 
phtude. This can be inferred from the analytical results 
in previous papers [l], 0, However, the strength of the 
squeezing effect on the two-identical particle correlation 



was not carefully investigated in those references. 

The analytical form of the HBT correlation function is 
obtained by substituting the chaotic amplitude 



Ge(ki,k2) 



(27r)i I 



|2d3 
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exp 



im{u)R 
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together with the spectrum given in Eq. (|4]) , into Eq. ^ . 
The finite emission time factor, multiplying the square 
modulus of Eq. (O, is now Fc{M) = [1 + (cji -a;2)^Af^]"\ 
For stressing the HBT effects in the (jxf) case, we selected 
the region (small Iq^jl) where the particle- antiparticle 
correlation is not significant and the HBT is relevant to 
Eq.©. 
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FIG. 3: The plots show the HBT correlation function in the 
absence of squeezing (top) and when it is present (bottom). 

This is seen in Fig. |31 The top part shows the ef- 
fect of radial flow alone on the HBT correlation function, 
while in the bottom, the joint effects of flow and squeez- 
ing are shown. We see that, without squeezing, the flow 
broadens the correlation curves, as expected, since the 
expansion reduces the size of the region accessible to in- 
terferometry. When the squeezing effects are present, 
they seem to oppose to the flow effects, almost canceling 
the broadening of the correlation function due to flow for 
large |K|, another striking indication of in- medium mass 
modification. 

CONCLUSIONS 

In this work we suggest an effective way to search for 
the squeezed bosonic correlations in heavy ion collisions 



at RHIC and, soon, at the LHC. We argue that the suit- 
able variable to experimentally search for the squeezed 
correlation function is the average momentum of the 
pair, 2|Ki2|, the non-relativistic limit of the relativis- 
tic variable, Qbbc = 2-^ {ijJiuj2 — Kt^K^^) [gj]. We show 
that, in the presence of flow, the signal is expected to 
be stronger over the momentum regions shown in the 
plots, i.e., roughly for < |2K| < 100 MeV/c (depend- 
ing on R) and 500 < |q| < 1500 - 2000 MeV/c, sug- 
gesting that flow may enhance the strength of the BBC 
signal, facilitating its experimental discovery. Another 
important result found within our simple non-relativistic 
model is that the squeezing could also distort the HBT 
correlation function, leading to effects opposing those of 
flow, almost neutralizing it for large values of jK^^]. For 
emphasizing the dramatic effects induced by in-medium 
hadronic mass modiflcation on the correlation functions, 
we chose a constant mass-shift that leads to the maximal 
intensity, based on results of Fig. ([2]). For cfxf) mesons, 
this corresponds to m* w 1 GeV, roughly a 2% reduc- 
tion in the (/) mass, as compared to its asymptotic mass 
(m ~ 1.02 GeV). As stressed before, a more realistic 
treatment should consider a detailed prescription for the 
mass modiflcation, based on models that predict its de- 
pendence on the particles' momenta and its distribution 
in the hot and dense system. 



The above procedure is also applicable to other parti- 
cles, such as kaons. The corresponding results are dis- 
cussed in Ref. [l0|. Finally, it is important to note that 
all the effects shown here should exist only if the particles 
have their mass modified in the hot and dense medium. If 
no modification happens, the squeezed correlation func- 
tions would be flat unity, and the HBT correlation func- 
tions would behave as usual. However, if the particles' 
masses are indeed modiflcd, the experimental discovery 
of squeezed particle-antiparticle correlation (and the dis- 
tortions pointed out in the HBT correlations) would be 
an unequivocal signature of in-medium modifications by 
means of hadronic probes! 
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